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Abstract 

We set a preliminary 95% C.L. exclusion on the oscillation frequency of — mixing 
using a sample of 400,000 hadronic decays collected by the SLD experiment at the 
SLC during the 1996-98 run. In this analysis, B^ mesons are partially reconstructed by 
combining a fully reconstructed Ds with other B^ decay tracks. The Dg decays are recon- 
structed via the (pir and K*K channels. The 6-hadron flavor at production is determined 
by exploiting the large forward-backward asymmetry of polarized bb decays as well 

as information from the hemisphere opposite to the reconstructed B decay. The flavor of 
the Bg at the decay vertex is determined by the charge of the Dg. A total of 361 candi- 
dates passed the final event selection cuts. This analysis excludes the following values of 
the B^ — B^ mixing oscillation frequency: Arris < 1-5 ps~^, 2.6 < Am^ < 4.9 ps~^, and 
10.8 < Arris < 13.5 ps~^ at the 95% confidence level. 
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1 Introduction 



The Standard Model allows ^ oscillations to occur via second order weak inter- 
actions. The frequency of oscillation is determined by the mass differences, Am, between 
the mass eigenstates in the B^ system. The mass difference in the B^ system (Am^) and 
in the B^ system (Am^) are proportional to the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements \Vts\ and \Vtd\, respectively. A measurement of Am^ can in principle be 
used to extract the CKM matrix element \Vtd\- However, the extraction of \Vtd\ from Am^ 
is complicated by a large theoretical uncertainty on the hadronic matrix elements. The 
complication can be circumvented by taking the ratio of Am^ and Am^. In the ratio, the 
theoretical uncertainties partially cancel, giving 

Arrid m^o 

d 

where is estimated from Lattice QCD calculation to be (l.llj;Q;o4)^ [I], ^ and m^o (m^o) 
is the mass of the B^ (5^) meson. Therefore, a direct measurement of Arris, combined 
with the current measurement of Arrid, can be translated to a value of \Vtd\ with a 5% 
precision. This has motivated many experiments to search for B^ oscillations. As of now, 
all analyses have failed to observe significant signal and only lower limits on Arris have 
been set. 

A mixing analysis requires several key ingredients: a B^ enriched event sample, 
knowledge of the B^ meson flavor at the production and decay vertex, and the B^ proper 
decay time. In this paper, an analysis to search for B^ oscillations using fully reconstructed 
Ds is presented. We begin in section 2 with a brief description of the experimental 
apparatus. In section 3, the details of the event selection and Ds reconstruction are 
discussed. The two Ds decay channels used in the analysis are: 

B''^^D;+X, D;^4>n-, ct>^K+K-; (2) 

5° ^ D; + X, D; K*°K-, K*^ K+n-. (3) 

The next two sections describe the method for reconstructing the B^ proper decay time 
and determining the flavor of the 5° at the production and the decay vertices. Finally, 
the fitting procedure and the results are discussed in the remaining sections. The re- 
sult presented in this paper is based on the data collected during the 1996-1998 runs 
which consist of 400,000 hadronic decays with an average longitudinal electron beam 
polarization of about 73%. 

2 The Detector 

The Stanford Large Detector (SLD) is a general purpose particle detector designed to 
study the decay of boson produced at the Stanford Linear Collider (SLC). A detailed 
description of the detector can be found here P, §]. In this section, the main features of 
the detector components that are important to the analysis are summarized. 
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At the heart of the SLD is a CCD pixel vertex detector (VXDS) with over 300 
miUion pixels. It consists of 3 concentric layers at radii of 2.8, 3.8 and 4.8 cm. The polar 
angle coverage of the vertex detector extends from cos{6) of -0.85 to 0.85. The single hit 
resolution is roughly Afim in each coordinate and the measured track impact parameter 
resolutions are: 

33 33 
ar4> = 7.8 © - — — TTT^/im , arz = 9.7 © - — -^nr^l^m , 

where z axis points along the beampipe and the track momentum [p) is measured in 
GeV/c. The VXD3 is surrounded by the Central Drift Chamber (CDC) that extends from 
the inner radius of 20cm to the outer radius of 100cm. Reconstructed tracks from the 
CDC are linked with tracks from the vertex detector to improve the momentum resolution. 
The momentum resolution for the combined tracks (CDC+VXD3) is measured to be 

(yp^/pL = 0.0095 © 0.0026px , 

where p±{GeV/ c) is the momentum transverse to the beampipe. The next layer of detector 
is the main particle identification system at SLD, the Cherenkov Ring Imaging Detector 
(CRID). The CRID comprises two radiator systems (C6-F14 liquid and C5F12 gas) for K/ n 
separation over a large momentum range (0.35-35.0 GeV/c). The barrel CRID provide 
particle identification over the central 70% of the solid angle. Immediately outside the 
CRID is the Liquid Argon Calorimeter (LAC), which is used primarily for triggering, en- 
ergy reconstruction, and electron identification. The calorimeter towers are made of planes 
of lead plates separated by non-conducting spacers and immersed in liquid argon. The en- 
ergy resolutions for electromagnetic showers is measured to be a^m./E = 15%/ ^E{GeV), 



and for the hadronic showers is estimated to be ahad/E = 60%/y i?(GeV^). 

The above four central sub-detectors are inside a solenoidal magnet that provides 
a uniform axial magnetic field of 0.6 Tesla. The outermost detector is the Warm Iron 
Calorimeter (WIC). The WIC consists of 14 layers of 5cm thick iron plates instrumented 
with streamer tubes between layers for muon identification. 



3 Event Selection and Ds Reconstruction 
3.1 hh Event Selection 

Events used in the analysis are first required to pass the hadronic event selection. An event 
is selected as a hadronic candidate if it satisfies the following conditions: has at least 7 
good CDC tracks (each track is required to have momentum transverse to the beampipe 
greater than 200MeV) that pass within 5 cm of the interaction point (IP), the thrust 
axis is within \cos{Othrust)\ < 0.85, and the total charged track energy is greater than 
18 GeV. The hadronic event selection removes essentially all di-lepton events {e^e" 

and other non-hadronic backgrounds. To further enhance bb events in the 

sample, hadronic candidates are required to have at least one topologically reconstructed 
secondary vertex ||^ with a vertex mass greater than 2 GeV that is corrected for the 
missing transverse momentum to partially account for the neutral particles. A neural net 
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is used in the vertex mass reconstruction to enhance the separation between bb and other 
qq events in the sample 0. Figure |l| shows the pt corrected vertex mass distribution 
for data and MC. A minimum vertex mass cut at 2 GeV yields an event sample with b 
hadron purity of about 98% and single hemisphere b tagging efficiency of about 56%. 

3.2 L>7 (t)T^-,K*^K- Reconstruction 

candidates are reconstructed by first pairing oppositely charged tracks to form a 
{K*^) candidate for the ^vr" {K*^K~) mode. A third track is then attached to form a Dg 
candidate. Charged tracks used in the combination are required to have at least 23 CDC 
hits (out of possible 80 hits), at least 2 vertex detector hits, and a combined x^/d.o.f. 
for the CDC+VXD3 fit of less than 8 to ensure good reconstruction. To maximize the 
discrimination between true Ds and combinatorial background events, kinematic infor- 
mation for the Ds candidate is fed into a neural net. The neural net inputs for the 07r 
[K*^K) mode include: K^K~ (/T+tt") invariant mass, fitted vertex probability of the Dg., 
total momentum of the D^, helicity angle 9* (angle between the {K*^) and the Dg fiight 
directions in the Dg rest frame), helicity angle A* (angle between the charged daughter of 
the Dg and the from {K*^) decay in the rest frame of the neutral meson), and par- 
ticle ID information for the three tracks. The complete list of neural net inputs is shown 
in Table |l|. The neural net is trained on Monte Carlo events generated using JETSET 7.4 



07r Mode 


K*^K Mode 


Dg vertex prob 


Dg vertex prob 


Pptot{Dg) 


Pptot{Ds) 


KK opening angle 


Kn (from K*^) opening angle 


Dg normalized decay length 


Dg normalized decay length 


helicity angle A* 


helicity angle A* 


helicity angle 9* 


helicity angle 9* 


particle ID of three tracks 


particle ID of three tracks 




Average Normalized 3-D impact 




parameter of KKtt tracks 



Table 1: Dg Neural net inputs - left column is for 07r and right column is for K*^K 

[H with full detector simulation based on GEANT 3.21 The neural net outputs are 
shown in Figure 0. 

The optimal neural net cut that maximizes the sensitivity of the analysis to -8° 
mixing is determined separately for each of the two Dg decay modes. The minimum cut 
for the 07r mode is determined to be at 0.9 and for the K*^K mode, is at 0.8. Furthermore, 
for the K*^K mode, the kaon from the K*^ decay is required to be identified by the CRID 
in order to suppress combinatorial and other non-D^ backgrounds. The KKtt invariant 
mass spectra are shown in Figures |^ and ^. The Dg mass peak is fitted separately for 
events with and without definite kaon ID and for Q=0 and Q=±l events, where Q is 
defined as the total charge of all tracks associated with the B decay. The details of the 



4 




Figure 1: pt corrected vertex mass distribution for data (red square) and MC (solid line). 
White histogram is bb events, blue histogram is cc events, and green histogram is uds 
events 




0.9 0.925 0.95 0.975 1 0.8 0.85 0.9 0.95 1 



NNOUT NNOUT 

Figure 2: Neural Net outputs for (a) (pre and (b) K*K^ modes. Dots are data, open 
histograms are Monte Carlo and solid black histograms represent MC simulation of true 
Ds events. 
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B vertex reconstruction and the rriKK-K niass fits are described in later sections. The 
estimated number of Dg candidates and combinatorial fractions are given in Table |^ and ^. 





Without De 
Q=0 


finite Kaon ID 
Q=±l 


With Defini 
Q=0 


te Kaon ID 
Q=±l 


^ of hadronic candidates 
# of semileptonic candidates 
Average combinatorial fraction 


101 

8 

42.3±2.1% 


57 
6 

53.8±2.7% 


54 
12 

21.7±2.2% 


40 
2 

14.1±1.4% 



Table 2: (pn mode- number of Ds candidates within ± 40MeV of the nominal mass and 
the estimated average combinatorial fractions. 





1 Kaon ID 
Q=0 


2 Kaons ID 
Q=0 


# of hadronic candidates 
jj^ of semileptonic candidates 
Combinatorial fraction 


40 
6 

32.9±3.3% 


30 

5 

22.3±2.2% 



Table 3: K*^K mode- number of Dg candidates within ± 40MeV of the nominal mass 
and the estimated average combinatorial fractions. 



4 B Proper Time and Vertex Charge Reconstruction 

The oscillation frequency is expected to be large in the Standard Model and therefore 
a time dependent study of the phenomenon requires precision reconstruction of the 5^ 
proper decay time. The proper decay time is calculated from the reconstructed decay 
length (/) and the momentum of the 5°: 

rB. = ^. (4) 

PIC 

where /?7, often referred to as "boost" , is the ratio between the momentum of the 5° and 
its mass. In the following sections, the algorithms for reconstructing the decay length and 
the boost are described. 



4.1 Decay Length Reconstruction 

The decay length of the is defined as the distance between the IP and the i?° decay 
vertex. The IP position is found by vertexing all tracks in the event to a common vertex. 
Taking advantage of the stability of the SLC beamspot, the IP position in the x-y plane 



6 




Figure 3: Distribution of KKn invariant mass for the (j)Ti mode. Plots in the left column 
(al,a2) are for the neutral candidates and plots in the right column (hl,h2) are for the 
charged candidates. The 07r sample is further subdivided into events with loose kaon ID 
(al,bl) and events with hard kaon ID (a2,b2). 




1.7 1.8 1.9 2 2.1 2.2 

mkkn (GeV) 




mkkn (GeV) 



Figure 4: Distribution of KKtt invariant mass for the K*K^ mode, (a) for events with 
1 kaon ID and (b) for events with 2 kaons ID. 
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(transverse to the beampipe) is averaged over 30 hadronic events. The position in z 
(beampipe direction) is calculated on an event-by-event basis. The estimated IP resolution 
for —>■ bb events is about Afim in x-y and 20/im in z. 

The B'^ decay vertex is found by vertexing the Dg track with other B decay tracks 
in the hemisphere. This is accomplished in two steps. The first step involves the selection 
of the seed vertex (preliminary estimate of the B^ decay vertex). To find the seed, the Dg 
track is individually vertexed with each quality track (excluding Dg daughters) in the same 
hemisphere, and the vertex that is farthest from the IP and upstream (or consistent with 
being upstream within 5a) of the Dg and has a vertex fit of Isss than 5 is chosen as the 
seed. Step two involves the separation of tracks into secondary decay and fragmentation 
tracks. The discriminating variable used is the L/D parameter, where L is defined as the 
distance from the IP to the point of the closest approach of the candidate track to the 
seed vertex axis(line joining IP to the seed), and D is the distance from the IP to the seed. 
A track is chosen as a secondary B decay track if it satisfies the following two conditions: 
1.) track L/D is greater than 0.5 and 2.) forms a good vertex with the Dg track (fit ^ 
5). The latter condition is imposed to reject spurious and B daughter tracks (from double 
charm decays) that do not point back to the B vertex. Finally, the selected tracks are then 
vertexed together with the Dg to obtain the best estimate of the B decay position. The 
resulting B decay length resolution is highly dependent on the decay topology. The list 
of decay length resolutions estimated from Monte Carlo for the various decay categories 
is shown in Table 



Decay Category 


Q=0 

Core aLifi'Tn) Tail aiifJ'm) 


Q=±l 

Core aLifi'Tn) Tail aiifJ^m) 


B^g DgX (right-sign) 


47 144 


51 184 


i?^ D^X (wrong-sign) 


89 292 


89 292 




70 271 


85 412 


Bu ^ D^X 


99 435 


72 236 


BBaryon — > DfX 


84 221 


84 221 



Table 4: Decay length resolutions for various decay topologies. Resolutions are param- 
eterized separately for neutral and charged B events except for wrong-sign Bg and B 
Baryon events. The resolution is parameterized by the sum of two gaussians with core 
fraction fixed to 60%. 



4.2 Boost Reconstruction 

To obtain the relativistic boost, the energy of the B meson has to be reconstructed. The 
total energy of the B meson is the sum of the energy of the charged and neutral daughter 
particles: 

^ — ^charged + ^neutral- (5) 

The charged energy is determined by summing all the charged tracks associated with the B 
decay assuming pion mass (except for the two kaons from the Dg decay). To estimate the 



8 



neutral energy contribution, five different techniques are used. The first four techniques 
are calorimetery based and use various constraints (beam energy, jet energy, mass 
and calorimeter information) to estimate the neutral energy of the B meson The 
fifth technique is based only on the kinematics of the decay (B vertex axis, charged track 
momentum and mass constraint) to estimate the neutral energy The results from 
the five algorithms are then averaged, taking correlations into account, to obtain the total 



B energy. The relative boost resolutions a 



are shown in Table ^ 





Q- 


=0 


Q= 


=±1 


Decay Category 


Core crfSj{%) 


Tail cTi3y{%) 


Core (T/3^(%) 


Tail (Tf3j{%) 


5° D-X (right-sign) 


7.9 


19.1 


10.1 


26.5 


B^ — >■ Dfx (wrong-sign) 


8.6 


19.7 


11.1 


28.3 


B'l ^ Dfx 


9.5 


19.6 


9.1 


24.7 


Bu ^ Dfx 


10.4 


30.0 


8.4 


21.1 


BBaryon — > DfX 


10.0 


25.4 


11.1 


39.6 



Table 5: Relative boost resolutions for various decay topologies. The resolution is param- 
eterized by the sum of two gaussians with core fraction fixed to 60%. 



4.3 Charge Reconstruction 

Nominally, the charge of the b hadron is the sum of the charge of the quality tracks 
associated with the B decay. To improve charge purity, tracks with only VXD3 hits (VX- 
alone vectors) that are not used in vertex and momentum reconstruction are also included 
in the B charge determination 0. For the (pn mode, both the neutral B candidates (Q=0) 
and the charged candidates (Q=±l) are used. The B^ purity in the charged sample is 
considerably lower than in the neutral sample and therefore has reduced weight in the fit. 
For the K*^K mode, only the neutral candidates are included in the final event sample. 
The reconstructed vertex charge distributions for data and Monte Carlo are shown in 
Figure |^. 



5 Flavor Tagging 

To decide whether mixing has occured for a given event, the flavor of the B^ has to be 
determined at the production and the decay point. In this section, the methods used to 
tag the initial and final state of the B^ are described. 

5.1 Final State Flavor Tag 

A significant fraction of the B^ decays result in a Dg in the final state. This knowledge 
is used not only to enhance the B^ purity in the data sample but also as a means to 
determine the flavor of the B^ at the decay vertex. For the final state tag, we assume the 
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Reconstructed B Vertex Charge 



Figure 5: Reconstructed B vertex charge distribution for data (dots), M.C. neutral B 
events (open histogram), and M.C. charged B events (solid histogram). The distribution 
is for the (pre and the K*^K modes combined. 

Ds to come from the b c transition which imphes that a B'^ would decay into a Df and 
a B^ would decay into a D~ . However, roughly 10% of the time the virtual W produces a 
Ds with an opposite sign. This process is the dominant source of final state mistag. The 
issue of final state mistag will be further addressed in a later section of this paper. 



5.2 Initial State Flavor Tag 

Several techniques are used to determine the initial state of the 5°. The most powerful 
method, unique to the SLD, is the polarization tag. In a polarized — > bb decay, the 
outgoing quark is produced preferentially along the direction opposite to the spin of the 
boson. Therefore by knowing the helicity of the electron beam and the direction of 
the jet, the flavor of the primary quark in the jet can be inferred. The analyzing power 
of the polarization tag is highly dependent on the polar angle of the jet w.r.t. the beam 
axis and the electron beam polarization. The probability that the jet is a b quark jet is 



P{b) 



l + A 



FB 



(6) 



where ApB '^^ the polarized forward-backward asymmetry, defined as. 



A 



FB 



. Ae 



cosdi 



1- A^P^I + cosWt 



(7) 



with Af) and A^ assumed to have the Standard Model values of 0.935 and 0.150, respec- 
tively. The polar angle d-r in equation (|^) is defined as the angle between the thrust axis, 
which points in the event hemisphere, and the electron beam direction. Pg is the electron 
beam polarization (Pg > for right handed and < for left handed electron beam). 
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In addition to the polarization tag, the momentum weighted jet charge technique 
is also used to determine the initial state of the B^. The quantity Qjet is defined as the 
sum of the charge of the tracks in the opposite hemisphere weighted by the longitudinal 
momentum of the track w.r.t. the thrust axis, namely 



trks 



Qjet = Y.\Pi-^ ' (8) 

i 

— * 

where T is the thrust axis and k (equal to 0.5) is a parameter determined from Monte 
Carlo that maximizes the jet charge separation between q and q jets. The probability of 
tagging the b quark hemisphere, given the opposite hemisphere jet charge Qjet, is 

P(6) = L^, (9) 

where the coefficient a is determined from Monte Carlo to be -0.27. 

To further enhance the initial state tag purity, additional sources of information 
from the opposite hemisphere are used. This includes: vertex charge, charge of lepton, 
charge of kaon and dipole charge [0]. All the available tags in a given event are combined 
to obtain the overall initial state tag probability. For this analysis, the average initial state 
correct tag probability is about 75%. 



6 Oscillation Studies 

The final step in a mixing analysis is to fit for the oscillation frequency. This entails 
a search for periodic oscillations in the proper decay time distribution of events tagged 
either as mixed or unmixed. We expect the proper time distribution, in the limit that the 
lifetime difference between the two mass eigenstates is small, to have the following time 
dependence for the mixed and unmixed events: 

Pmixed{r) = — (1 - cos{AmsT)) , (10) 

Zt 

Punmixed{r) = — (1 + cos(Am^r)) . (11) 

Zt 

Where tb^ is the lifetime and Arris is the mass difference, as introduced in the earlier 
section. The fitting function used in the analysis also includes the effect of detector 
resolution, reconstruction efficiency, mistag, and background events. These details are 
addressed in the next few sections. 

6.1 Amplitude Fit 

The amplitude fit [jl3| is a fitting method that is used in this analysis. The method, 
in essence, transforms the traditional likelihood fit into a "Fourier-like" analysis. In an 
amplitude fit, the likelihood function is modified by introducing a term A, the amplitude, 
in front of the cosine terms (cos(Amsr) Acos(Amsr)). Instead of fitting for Am,; 
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directly, Am^ is fixed to a particular value and the parameter A is fitted for. The 
fit for A is repeated for a range of Am^ values to produce the amplitude plot. The 
amplitude is expected to be consistent with zero for values of Am^ far from Am*^"^ and 
to reach unity at the true mass difference value. If the oscillation frequency is large and 
no signal is observed, the range of Arris for which A + 1.645o"^ < 1 can be excluded at the 
95% confidence level. The 95% C.L. sensitivity is defined as the value of ArUg at which 
1.645aA = 1. 

6.2 Likelihood Function 

To perform the amplitude fit, we first need to construct a likelihood function that describes 
the proper time distribution of the events in the data sample. The events in the final 
sample can be divided into seven main sources, each with its own proper time distribution 
function, F^. The seven physics sources are: 

• F^^' = 5^ — > D~X right sign decays + c.c. 

• Fq^- = Bg wrong sign decays + c.c. 
. Fb, = ^ Dfx + c.c. 

• Fb± = S± ^ Dfx + c.c. 

• Fbb = B Baryon DfX + c.c. 

• Fee — primary charm quark — > DfX + c.c. 

• Fcornb = Combinatorial events. 

The proper time distribution of the event sample is the sum of the seven physics functions 
with the contribution from each source weighted by its fraction in the sample and the 
appropriate normalization constant. The resulting normalized probability distribution 
function for the mixed events is 

Cjr.s. jw.s. 

^-^Fbb + ^Fe)i + [1 - fDsiniKK.)] ■ Feomb , (12) 

where Tree is the reconstructed proper time, fn^ is the fraction of Ds in the sample, fx is 
the fraction of category x in the Dg signal peak, and Ni is the normalization constant for 
category i. The physics functions for the mixed events have the following form: 

FBsi^ree) = / [{1 - Vi) - ^ ■ COs{AmsT)) + Tji {1 + A ■ COs{AmsT))] ■ 

ei(r) • Gi(r,r,.ec) dr. (13) 
FsfiTree) = / -7, [(1 - 7?^) (1 + A • cos(Am,T)) + (1 - ^ • cos(Am,r))] • 
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e2{r) ■ G2{T,Trec) dr. (14) 

POO g-'^/'^Sd 

FB^Trec) = / — [(1 - T/fiJ (1 - cos(Amdr)) + r^fi, (1 + cos(Amrfr))] ■ 

Jo 2tb^ 

^sir) ■ G3{r,Trec) dr. (15) 

FB±iTrec) = / VB^ ' (^4.{t) ' Gi^T, Tree) dr. (16) 



OO p-'^hsf 



FBB{Tree) = / ^BB " £5(7") ' G'5(r, T^ec) C??-. (17) 

Jo ZTBbb 

Fee{Tree) = 5.47e-^'-'=^/°-2^°^ * er/( 18.35w) (parameterized /rom M.a). (18) 

Fcomb{Tree) = {distribution taken from data sidebands). (19) 

Where, r is the true proper time, rji is the initial state mistag probabihty, rjB^^g are 
the overall mistag probabilities for Bd,±^B events, and the last two functions: ej(r) and 
G'i(r, Tree) are the vertex efficiency and the resolution functions for category z, respectively. 
The proper time distribution for the combinatorial events is taken directly from the data 
using events in the sidebands. Detailed M.C studies have shown that the combinatorial 
distribution is well modelled by the sidebands and that using the sideband distribution 
does not introduce a bias in the amplitude fit. The corresponding probability distribution 
for the unmixed events is obtained by replacing rj with (1 — 77) in the physics functions. 
Finally, the unbinned likelihood function is defined as the product of the probabilitiy of 
the mixed and unmixed events 

mixed unmixed 

L W Pmxiedij'ree) W Punmxiedi^Tree) ■ (^0) 
i 3 

6.2.1 Data Compositions and Event Mistag Rate 

The Dg fraction (/dJ in equation (|12|) is the fraction of true Dg events in the sample. 
Instead of using the average value from the mass plot, the fraction is determined on 
an event-by-event basis using the reconstructed Dg mass of the candidate event. Events 
close to the nominal Ds mass are more likely to be true Dg events than combinatorial 
events, therefore have more significance in the analysis. The Dg signal and background 
parameterizations used in the /d^ calculation are taken directly from the rnxKn mass 
plot. The fitting function for the Dg mass peak is the sum of two gaussians with the 
same mean and 60% core fraction. The widths used in the individual fits are fixed to 
the widths from the combined fit. The background is parameterized by a second order 
polynominal (combinatorial background) and a gaussian function(Z}"'" mass peak). For 



the gaussian, only the amplitude is allowed to float; the mean is fixed to 99.2MeV [|14 
below the fitted Dg peak and the width is taken from the Monte Carlo distribution. The 
six dominant sources of Dg production are outlined in the beginning of this section. The 
relative fractions of the sources are determined from the Monte Carlo and are estimated 
separately for the charged and neutral samples as well as for hadronic and semileptonic 
decays. The semileptonic decays are defined as events where a B decay track (not from 
the Dg decay) is identified as either a muon or an electron. The relative fractions used in 
the fit are listed in Table || and ^ 
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Hadronic Decays 


Semileptonic Decays 




Q=0 


Q=±l 


Q=0 


Q=±l 


-f r.s. 


0.556 


0.283 


0.822 


0.452 


fw.s. 

JB, 


0.066 


0.034 


0.036 


0.020 


fSd 


0.260 


0.160 


0.122 


0.129 


fB± 


0.046 


0.452 


0.013 


0.387 


Ibb 


0.053 


0.052 


0.007 


0.011 


fee 


0.019 


0.019 


0. 


0. 



Table 6: Fractions for 07r mode. 





Hadronic Decays 
Q=0 


Semileptonic Decays 
Q=0 


fr.s. 
JBs 


0.558 


0.766 


fw.s. 

JBs 


0.055 


0.075 


fSd 


0.252 


0.117 




0.048 


0.014 


Ibb 


0.051 


0.028 


fee 


0.036 


0. 



Table 7: Fractions for K*^K mode. 



The contribution to the likelihood function is divided into the right-sign and 
wrong-sign decay terms, and therefore, by construction, the final state mistag rate is 
and 1, respectively, and does not enter explicitly into the likelihood. Instead, the effective 
final state mistag is accounted for in the relative 5° right-sign and wrong-sign fraction. 
For the contribution, the right-sign and wrong-sign decays are not treated separately 
and the overall event mistag rate factors in both the inital state {r]i) and final state [f]^'^) 
mistags tibj^ = rii{l — rj^'') + f]f'^{^ ~ Vi)- ^ind B Baryons, the overall event mistag 

rate is obtained by subsituting rj^'' in the previous equation with r]^^ and t]^^ . The final 
state mistag rates are listed in Table M and the various measured branching ratios used 



in the calculation are listed in Table 10 



6.2.2 Vertex Efficiency, Resolution Functions and Normalizations 

To complete the discussion on the likelihood function, two remaining effects have to be 
addressed. First, the vertex efficiency function is included in the physics functions to 
account for the lower efficiency for reconstructing events with short proper time (events 
close to the IP). The efficiency functions are taken from the Monte Carlo simulations and 
are parameterized by the function: 

<-) = ^i^^ + ^3. (21) 
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rjf (hadronic decay) 


rjf (semileptonic decays) 




79.9% 


88.5% 




77.4% 


93.8% 


B Baryons 


92.1% 


88.9% 



Table 8: Final state mistag rates for the hadronic and the semileptonic (events with a 
lepton attached to the B vertex) decays. For the B^, B^, and B baryon events, the 
primarily comes from the W decay, therefore, by definition, the final state mistag rates 
for those events are greater than 50%. 



Equation (^) models the vertex efficiency well for the 5°, -B^ and B Baryons events. 
However, for the 5^, the vertex efficiency is also expected to decline at large proper 
times due to the B vertex charge requirement on the final sample (only candidates with 
Q=0,±1 are kept). Therefore an additional exponential term {P^e'^'''^) is added to the 
function to model the behavior of the B^ events. The coefficients (-P1-5) are listed in 
Table | 



Category 


Pi 


P2 


P3 


P4 




ei,e2 


0.148 


-6.856 


0.033 






€3 


0.037 


-4.448 


0.032 






£4 


0.061 


-0.985 


-1.209 


1.243 


0.007 


es 


0.056 


-3.532 


0.002 







Table 9: Parameterizations for the vertex efficiency function. 



The second effect concerns the proper time resolution. The effect is accounted 
for by introducing a resolution function in the convolution integrals, as shown in equa- 
tions (|l^) to (0). The proper time resolution can be expressed in terms of the decay 
length and boost resolutions ((Xi,,^^/?) as in 




1/2 



(22) 



where the indices (i and j) refer to core or tail. The proper time resolution contains a 
constant term that depends on ai and a term that rises linearly with proper time that 
depends on cr^^. This behavior is illustrated in Figure ^ for the B^ events. In practice, 
there are four ar distributions given by the various and cx^^ core-tail combinations. 
The resolution function is obtained by summing all four contributions 



G(r, Tree) = f!G^ + mC^t + Gt,) + f.'Gu. 



(23) 
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'0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
Proper Time (psec) 



Figure 6: a^- distributions for D~X events. The lower red line is calculated using 

(Ti(core) and a.yp{core). The upper blue line is calculated using aiitail) and a^pitail). 
The background histogram is the proper time distribution of the data events. 



with /c=0.6, /t=0.4, and Gjj defined as 



The complete likelihood function is normalized by dividing each physics term by 
a normalization constant. The normalization constant (A^j) is calculated for category i 
by integrating the sum of the mixed and unmixed physics functions for source i over all 
reconstructed proper times, 

/ + 00 
^prnixed ^ punmixed^ (25) 
-oo 

The normalization is required to ensure that the weight of each physics source in the 
likelihood function is not biased by the efficiency and the proper time resolution of the 
event . 



6.3 Ds+Tracks Amplitude Fit Results 



The physics parameters used in the amplitude fit are listed in Table |T0|. The systematic 
uncertainties are calculated based on the formula 

C^new ^nominal \ 
A ^ A I f26) 

^nominal ] ' ^ ' 

(^A / 

The physics parameters are varied by ±lcr in the fit to obtain the systematic uncertainties. 
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Parameter 


Value and uncertainty 


Reference 


^BBaryons 


1.464 ± 0.057 
1.562 ±0.029 ps 
1.656 ± 0.025 ps 
1.208 ±0.051 ps 
0.476 ±0.016 ps-^ 




15 
15 
15 
15 
15 


1 
1 

1 

1 


fib - fi,") 
f{b — > B Baryon) 


0.1000 ±0.012 
0.4010 ±0.010 
0.0990 ±0.017 




15 
15 
15 


1 
1 

1 

1 


R, ■ Bib ^ 5°) • i3(i?o ^ D+X) . ^ 07r+) 
S(6 ^ ly- ^ D,-) ■ B{D; (/.TT-) 
^ Z}±X) ■ i3(Z}7 ^ 4>n-) 
B{Ba.u ^ Dfx) 
B(c ^ D7) ■ B{D; ^ 07r-) 


(6.2lllj:^^) X 10-^ 
(3.66 ±0.45) X 10-3 
(3.71 ±0.28) X 10-3 
0.172 ±0.083 

(3.4 ±0.3) X 10-3 




16 
15 
15 
15 


1 



Table 10: B lifetimes, Arrid, B production fractions, and various branching ratios assumed 
in the amplitude fit. The uncertainties for the branching ratios do not include uncertainty 
from Bt{Ds —>■ (pn). The right-sign Ds production fraction [Rh ■ B{b —>■ B^) ■ B{B'^ 
D^X) ■ B{Df — > (f)TT^)) is obtained by combining the direct measurement |jl6| with the 
semileptonic measurement W^ assuming factorization. 



In addition, the uncertainity on the Dg signal fraction (/dJ is estimated from the Ds mass 
fit and is varied by roughly 8%. The initial state tag probability is varied by ±0.02. The 
decay length resolution and the boost resolution uncertainties have not yet been studied 
in great detail and conservative estimates of 10% on ai and 30% on are used. A 
summary of the systematic uncertainties is given in Table |Tl] 

The resulting amplitude plot is shown in Figure |^ The measured values of the 
amplitudes are consistent with zero within the range of Arris considered and no evidence 
of a signal is observed. This analysis excludes the following values of the B^ — 5° mixing 
oscillation frequency: Arris < 1-5 ps"^, 2.6 < Am,, < 4.9 ps"^, and 10.8 < Arris < 13.5 
ps"^ at the 95% confidence level (C.L.). The sensitivity at the 95% C.L. is 1.4 ps^^. A 
comparison of the amplitudes and errors at Arris = 15 ps~^ for the various B^ mixing 
analyses is shown in Figure |^ (this analysis is listed as "SLD Ds" in the figure). 
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n- 
E 

< 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

SLD Ds+TRKS 

j data ± 1 a 
1.645 a -e- sensitivity 1.4 ps"^ 



I I data ± 1.645 a 

I I data ± 1.645 o (stat only) 



I ' ' ' ' I 




Am^ (ps" ) 

Figure 7: Ds+Tracks amplitude fit plot. 



ALEPH D Ji 
(91-9S) 

ALEPH B +D I 
(91-95, prelf) 

ALEPH 1 
(91-95, no D^l, prel.) 

CDF ld)/l 
(92-95) 

DELPHI B +D h 
(92-95, pret.) 

DELPHI D l+M 
(92-95, prel.) 

DELPHI 1 
(92-95, prel) 

DELPHI vtx 
(94-95, prel) 

OPALl 
(91-95) 

OPAL D I 
(91-9S) 

SLD 1+D 
(96-98, prel.) 

SLD dipole 
(96-98, prel.) 

SLD D, 
(96-98, prel.) 

World average (prel.) 



M 

'i i' 

-6 -4 -2 



amplitude (sensitivity) 
4.65 ± 3.74 1 1.07 ( 4.1 ps"') 


0.40 ± 1.31 ± 0.26 


(7.4ps-') 


1.49 ± 0.88 ± 0.31 


(11.7 ps-') 


-0.14 ± 2.00 ± 0.51 


( 5.1 ps-') 


0.45 ± 3.58 ± 1.93 


( 3.2 ps-') 


-0.54 ±1.90 ±0.32 


( 8.2 ps-') 


-1.59 ± 1.87 ± 0.19 


(7.8ps-') 


1.27 ± 3.27 ± 0.27 


( 5.0 ps-') 


-1.25 ±2.34 ±1.91 


( 7.2 ps-') 


-3.63 ± 3.05 t"4°2 


( 4.2 ps-') 


0.67 ± 1.07 


( 6.3 ps-') 


-1.20 ±1.09 to.!', 


(6.9ps-') 


0.78 ± 1.59 


(1.4ps-') 


0.25 ± 0.46 


(18.0 ps') 



B Oscillations 
Working Group 



amplitude at Am^ = 15.0 ps"' 

Figure 8: Amplitudes and errors at Ams=15ps~^ . 
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7 SLD Combination and Conclusion 



The analysis presented in this paper is complementary to the two more inclusive analyses 
at the SLD: Charge Dipole and Lepton+D. Charge Dipole and Lepton+D have a com- 
bined sensitivity of 10 ps~^ [0. Figure ^ shows the amplitude fit for all three analyses 
combined. The combined plot takes into account the correlated systematic uncertainties. 
Furthermore, the samples were selected such as to remove any statistical overlap between 
analyses. No evidence of a signal is observed up to Am^ of 25 ps~^ and the excluded 
regions at the 95% C.L. are: Am^ < 7.6 ps~^ and 11.8 < Am^ < 14.8 ps~^. The SLD 
combined sensitivity at the 95% C.L. is 13.0 ps~^. 
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Arris 

Measured amplitude A 
Statistical uncertainty (o"^*"*) 
Total systematic uncertainty (cr^^*) 



5 ps ^ 
-0.260 

±0.813 

+0.163 
-0.136 



10 ps-^ 
-0.106 

±1.058 

+0.129 
-0.098 



15 ps ^ 
0.781 

±1.591 

+0.161 
-0.193 



+0.U18 
-0.019 

-0.0006 
+0.0006 
-0.0007 
+0.0007 
-0.00003 
+0.00001 
-0.0007 
+0.0007 

+0.020 

-0.019 

+0.007 
-0.007 

-0.048 
+0.069 



+U.0U9 
-0.009 

-0.0005 
+0.0005 
-0.0002 
+0.0002 
-0.00002 
+0.00003 
-0.0002 
+0.0002 

+0.010 

-0.010 
+0.002 
-0.002 

-0.017 
+0.025 



+U.012 

-0.013 

-0.001 
+0.001 
+0.00002 
-0.00002 

+0.0001 
-0.0001 
+0.00008 
-0.00008 

+0.003 
-0.003 
+0.002 
-0.002 

-0.0006 
-0.0014 



20 ps-1 

-2.425 

±1.817 

+0.169 
-0.221 



+U.UU6 
-0.007 

-0.0002 
-0.0002 
-0.001 
+0.001 
+0.002 
-0.002 

-0.00008 
+0.00008 

-0.0007 

+0.0003 

-0.005 
+0.005 

-0.0068 
+0.0007 



'^B baryons 

Arud 

fib - 5°) 
f{b B Baryon) 

B(M DtX) . B{Dt - 
B{h ^W-^ D-) ■ B{D- 
B{Bd,u ^ Dfx) ■ B{D- - 
B{Bd,u ^ D-X)/B{Bd,u - 
B{c ^ D-) . B{D- 

Decay length resolution 
Boost resolution 

Initial state tag 



> (pTT ' 

^ 07r-) 
- DfX) 



+0.020 


+0.011 


+0.003 


-0.0008 


-0.020 


-0.010 


-0.003 


+0.0003 


+0.021 


+0.006 


-0.0009 


+0.008 


-0.021 


-0.006 


+0.0008 


-0.009 


+0.014 


+0.024 


+0.020 


+0.033 


-0.013 


-0.024 


-0.021 


-0.040 


+0.006 


+0.004 


-0.0004 


-0.005 


-0.006 


-0.004 


+0.0003 


+0.005 


+0.004 


+0.039 


+0.015 


+0.027 


-0.003 


-0.038 


-0.019 


-0.040 


+0.052 


+0.075 


-0.166 


-0.147 


-0.027 


-0.034 


+0.138 


+0.019 


+0.057 


+0.013 


+0.047 


+0.004 


-0.055 


-0.011 


-0.051 


-0.006 


-0.102 


-0.075 


-0.079 


-0.155 


+0.118 


+0.089 


+0.061 


+0.162 



Table 11: Table of statistical and systematic uncertainties for several Arris values. 
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Figure 9: Combined SLD amplitude plot. 
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